\ ^ \ TKHR8200M080 

APPLICATION OF NEURO-OCULAR WAVEFRONT DATA 
IN VISION CORRECTION 

STATEMENT RELATED TO COLOR DRAWINGS 

[0001] The patent or application file contains at least one drawing executed in color. Copies of 
this patent or patent application publication with color drawings(s) will be provided by the Office 
upon request and payment of the necessary fee. 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0002] This application claims the benefit of U.S. provisional patent application serial number 
60/458,480, filed on March 28, 2003, and having the title "Systems and Methods for Measuring, 
Correcting, and Optimizing Vision," which is incorporated herein by reference in its entirety. 

COPYRIGHT NOTICE 

[0003] A portion of the disclosure of this patent docimient contains material, which is subject to 
copyright protection. The copyright owner has no objection to the facsimile reproduction by 
anyone of the patent document or patent disclosure as it appears in the U.S. Patent and 
Trademark Office patent file or records, but otherwise reserves all copyright rights whatsoever. 

TECHNICAL FIELD 

[0004] The present disclosure relates generally to metrics and, more particularly, to visual 

metrics associated with vision correction. 
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BACKGROUND 

[0005] Traditionally, visual acuity has been measured using eye charts, in which subjects 

interactively provided vision-related information by reading characters of various sizes from the 
eye chart. The refractive error of the eye has traditionally been determined by having the subject 
recursively view objects through various sphero-cylindrical lenses, and selecting the combination 
of lenses that provides the clearest image. While such acuity and refraction measurements are 
still widely used, they provide only gross measurements of the visual system. 

[0006] Unlike, refraction, ocular aberrometry provides greater details with reference to the 
refractive properties of the eye. In ocular aberrometry, the refractive properties of the eye are 
measured using various lenses and projections, typically sampling multiple points across the 
pupil rather than a single measurement as performed during refraction, thereby providing more 
detailed information on higher-order optical anomalies. The aberration information provided by 
ocular aberrometers is designated herein as "ocular wavefront error." Examples of ocular 
aberrometers and ocular aberrometry are disclosed in U.S. Patent Numbers 5,777,719, 6,095,651, 
and 6,51 1,180, which are incorporated herein by reference in their entireties. 

[0007] While ocular aberrometers provide detailed information on the optical characteristics of 
the eye, these devices usually provide no information with reference to the neurological 
pathways between the eye and the brain and they are not capable of locaHzing the retinal plane 
preferred by the patient. In other words, the integrated visual pathway between the eye and the 
brain can introduce another transfer fiinction, which is undetectable by ocular aberrometry. 

[0008] In order to address these deficiencies, a technique known as visual aberrometry has been 
developed, in which the patient's preferred retinal plane is localized and the effects of the 
neurological pathways are taken into account by combining subject feedback with the physio- 
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optical characteristics of the eye. The aggregate effect of the neurological pathway in 
combination with the characteristics of the eye is referred to herein as "neuro-ocular wavefront 
error" or "neuro-ocular wavefront data." Examples of visual aberrometers are described in 
greater detail in U.S. Patent Numbers 5,258,791, 6,000,800, and 6,099,125, which are 
incorporated herein by reference in their entireties. 
[0009] Visual aberrometry, and the acquisition of neuro-ocular wavefront data, is still relatively 
new and continually improving. Along with the continued development of visual aberrometry, 
there is corresponding effort to exploit the full potential of visual aberrometry. 

SUMMARY 

The present disclosure provides systems and methods for diagnosing and treating subjects using 
neuro-ocular wavefront data. As such, in some embodiments, among others, neuro-ocular 
wavefront data is obtained, and one or more characteristics of a visual system are ascertained 
from the neuro-ocular wavefront data. 
[001 0] Other systems, devices, methods, features, and advantages will be or become apparent to 
one with skill in the art upon examination of the following drawings and detailed description. It 
is intended that all such additional systems, methods, features, and advantages be included within 
this description, be within the scope of the present invention, and be protected by the 
accompanying claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[001 1] Many aspects of the disclosure can be better understood with reference to the following 
drawings. The components in the drawings are not necessarily to scale, emphasis instead being 
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placed upon clearly illustrating the principles of the present disclosure. Moreover, in the 

drawings, like reference numerals designate corresponding parts throughout the several views. 
[0012] FIG. 1 is a diagram illustrating an eye, and various axes defined in relation to the eye. 

[0013] FIG. 2 is a diagram illustrating an embodiment of a pupillometry image, which may be 

used in a refi-actometer adapted to acquire neuro-ocular wavefiront data. 
[0014] FIG. 3 is a diagram illustrating an embodiment of a sampling matrix in a refractometer 

adapted to acquire neuro-ocular wavefi-ont data. 
[0015] FIG. 4 is a diagram showing, more clearly, the treatment envelope of FIG. 3. 

[0016] FIG. 5 is a diagram showing embodiments of various pupil-sampUng maps. 

[0017] FIG. 6 is a diagram showing an example of an operator's display console of a 

refi-actometer adapted to acquire neuro-ocular wavefi'ont data. 
[0018] FIG. 7 is a diagram showing another example of the operator's display console of a 

refi-actometer adapted to acquire neuro-ocular wavefi-ont data. 
[0019] FIG. 8 is a diagram showing a display console with a sample alignment target and a test 

spot. 

[0020] FIG. 9 is a diagram showing a coordinate system of the neuro-ocular wavefront data and 
a coordinate system for an ablation map that corresponds to the neuro-ocular wavefi-ont data. 

[0021] FIG. 10 is a diagram showing an optical zone and a transition zone associated with the 

ablation map of FIG. 9. 

[0022] FIG. 1 1 is a diagram showing an embodiment of a sampling map that has a predefined 

matrix of sampling points. 
[0023] FIG. 12 is a diagram showing an embodiment of a sampling map that permits 

customization of sampling points. 
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[0024] FIG. 13 is a diagram showing an embodiment of a subject data entry screen for a 

refiractometer adapted to acquire neuro-ocular wavefront data. 
[0025] FIG. 14 is a diagram showing an embodiment of a video display that facilitates on-screen 

measurements of physical characteristics of the eye. 
[0026] FIG. 1 5 is a diagram showing an embodiment of a video display that permits annotation 

of various features of the eye. 
[0027] FIG. 16 is a block diagram showing an embodiment of a system configured to aggregate 

data fi-om multiple refractometers. 
[0028] FIG. 17 is a flowchart showing an embodiment of a method for acquiring and processing 

neuro-ocular wavefront data. 
[0029] FIG. 18 is a flowchart showing, in greater detail, an embodiment of the data-acquisition 

step from FIG. 17. 

[0030] FIG. 19 is a flowchart showing, in greater detail, an embodiment of the step of 

identifying visual zones, from FIG. 18. 

[0031] FIGS. 20A and 20B are flowcharts showing, in greater detail, an embodiment of the step 

of interactively obtaining information, from FIG. 1 8. 

[0032] FIG. 21 is a flowchart showing, in greater detail, an embodiment of the step of correlating 

neuro-ocular wavefront data to vision parameters, from FIG. 17. 

[0033] FIGS. 22A and 22B are flowcharts showing, in greater detail, an embodiment of the data- 

analysis step from FIG. 17. 

[0034] FIG. 23 is a flowchart showing, in greater detail, another embodiment of the data-analysis 

step from FIG. 17. 

[0035] FIG. 24 is a flowchart showing another embodiment of a method for acquiring and 
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processing neuro-ocular wavefront data. 
[0036] FIG. 25 is a flowchart showing an embodiment of a method for acquiring and processing 

nexuro-ocular wavefront data, in which the processing of subsequently-acquired data is refined. 
[0037] FIG. 26 is a flowchart showing an embodiment of a method for acquiring and processing 

neuro-ocular wavefront data, in which multiple sets of data are statistically analyzed. 
[0038] FIGS. 27A and 27B are flowcharts showing another embodiment of a method for 

acquiring and processing neuro-ocular wavefront data, in which multiple sets of data are 

statistically analyzed. 

[0039] FIG. 28 is a table showing definitions of various variables related to the locations of 

sampling points and their respective corrections. 
[0040] FIG. 29 is a table showing definitions of the coordinate systems in FIG. 9. 
[0041] FIG. 30 is a table showing variable definitions used for computation of an ablation profile 

from acquired neuro-ocular wavefront data, as described with reference to FIG. 10. 
[0042] FIG. 3 1 is a table showing an example list of Zemike fimction terms. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

[0043] Reference is now made in detail to the description of the embodiments as illustrated in 

the drawings. While several embodiments are described in connection with these drawings, 
there is no intent to limit the invention to the embodiment or embodiments disclosed herein. On 
the contrary, the intent is to cover all alternatives, modifications, and equivalents. 

[0044] As noted above, traditional approaches to vision correction include measuring vision with 
visual acuity determined by eye charts and measuring the refractive error of the eye with test 
lenses. More advanced approaches employ ocular aberrometry, which provides greater details 
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on the refractive properties of the eye. While both of these approaches provide useful results for 
vision correction, they are nonetheless deficient in many ways. For example, visual acuity 
provides nothing more than a gross estimate of the performance of the patient's visual system. 
Refraction provides only the best sphero-cylindrical correction that is averaged across the pupil. 
Ocular aberrometry, while providing more detailed refractive characteristics of the eye, provides 
no information on the preferred retinal plane and no information related to the neurological 
pathways associated v^dth the eyes. 

[0045] An improvement to both refraction and ocular aberrometry is visual aberrometry, in 
which the patient preferred retinal plane and the effects of the neurological pathways are taken 
into account by incorporating a subject feedback in a testing process that assesses multiple 
points within the pupil. The aggregate effect of the location of the patient's preferred retinal 
plane and the effects of the neurological pathways, is known as a "neuro-ocular wavefront error" 
or "neuro-ocular wavefront data." The present disclosure provides various systems and methods 
for exploiting the full potential of visual aberrometry. In other words, the present disclosure 
provides systems and methods in which neuro-ocular wavefront data can be used for diagnosis, 
treatment, or a combination of both. 

[0046] In that regard, in some embodiments, neuro-ocular wavefront data is obtained from a 

subject (or, more specifically, a patient) and, thereafter, processed. In one embodiment, among 
others, processing of the neuro-ocular wavefront data includes correlating the neuro-ocular 
wavefront data with various subject^dependent or procedure-dependent parameters (hereinafter 
referred to in the aggregate as "vision parameters"). For example, the neuro-ocular wavefront 
data is often represented as an equation. Specifically, for vision correction, an example equation 
may be a Zemike polynomial with multiple coefficients. The coefficients in the Zemike 
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polynomial, either alone or in combination with other coefficients, represents a particular 
characteristic of the patient's visual system. 
[0047] In other embodiments, the correlated information is used for designing treatments. In 
some embodiments, the treatment can be designed by inverting the neuro-ocular wavefront data. 
For example, if the neuro-ocular wavefront error is defined by a Zemike polynomial, then the 
inversion of the polynomial provides one or more correction factors that can be applied, thereby 
effectively correcting for the apparent error that is introduced by the aberrations in the neuro- 
ocular system. For some embodiments, the inversion of the equation can be computationally 
derived using a least-squares fitting algorithm, a minimization algorithm, or other algorithms 
designed for curve fitting. Specifically, for neuro-ocular wavefront data that is represented by a 
Zemike polynomial, various coefficients of the Zemike polynomial may provide information on 
the near-vision or far-vision refraction sphere, refraction cylinder, and refraction axis. The 
curve-fitting algorithm provides information on how much correction will be needed to remedy 
the errors introduced by the refraction sphere, refraction cylinder, and/or refraction axis and the 
higher order aberrations of the visual system that are not corrected by sphero-cylindrical 
corrections. 

[0048] Having described general concepts associated with various embodiments of the 

invention, attention is turned to FIGS. 1 through 27B, which provide greater details related to 
various embodiments of systems and methods for acquiring and processing neuro-ocular 
wavefront data. 

[0049] FIG. 1 is a diagram illustrating an eye 1 10a, and various axes defined in relation to the 

eye 1 10a. Specifically, FIG. 1 shows the right eye 1 10b, the left eye 1 10c, and the nose 120 
from the viewpoint of the clinician (or other individual performing the test on the subject). The 
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coordinate system is right-handed insofar as the horizontal axis x extends from the left of the 
clinician (right of the subject) to the right of the clinician (left of the subject). The vertical axis y 
extends inferior to superior. Thus, the negative y-axis extends downward (toward the feet) while 
the positive y-axis extends upward (toward the head of the subject). 

[0050] The x-axis and the y-axis define the pupillary plane. Since the corneal plane is parallel to 
the pupillary plane, the x-y plane also defines the corneal plane. The z-axis is defined as the axis 
normal to the pupillary plane. In that regard, the direction of light propagation in the eye 1 10a is 
defined by the positive z-axis. Due to the spherical nature of the eyes 1 10b, 1 10c (hereinafter 
simply referred to as 1 10), the coordinate system is defined using polar coordinates, where r 
represents the radius from the x-y origin, and 9 represents the angle measured between the ray 
projected from (0, 0) to (x, y) with positive angles corresponding to counter-clockwise rotations. 
The normahzed radius, is computed by taking the radius, r, and dividing by a normalization 
radius, rnorm. In this document, in most cases, rnorm is taken to be 3.5 nmi. 

[0051] Additionally, FIG. 1 shows a compensation angle dy along the y-axis. This angle is the 
angular deviation that corresponds to a degree of refractive error, along the y-direction, for a 
given point on the eye. It should be appreciated that there is a similar compensation angle dx 
along the x-direction. The two compensation angles dx, dy provide an indication of a correction, 
as described in greater detail below. For simplicity, a summary of these variables is provided in 
FIG. 28. 

[0052] FIG. 2 is a diagram illustrating an embodiment of a pupillometry image, which may be 
used in a refractometer adapted to acquire neuro-ocular wavefront data. An example of an 
aberrometer adapted to acquire neuro-ocular wavefront data can be foimd in U.S. patent number 
5,258,791 (hereinafter, "the 791 patent"), having the title "Spatially Resolved Objective 
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Autorefractometer," filed on July 24, 1990, and issued on November 2, 1993, which is 
incorporated herein by reference, as if set forth in its entirety. For simplicity, the 
autorefi"actometer of the 791 patent is also referred to as the "InterWave" refi-actometer. 
[0053] As shown in FIG. 2, the pupillometry image includes an image of the eye 110, which 

includes a pupil 210. For proper orientation, both the x-axis and the y-axis are shown in FIG. 2. 
In addition to the image of the eye 1 10 and the pupil 210, the pupillometry image includes a 
pupil tracking box 220, which dynamically tracks the movement of the pupil 210. Thus, the 
motion of the pupil tracking box 220 is substantially synchronous with the motion of the pupil 
210. 

[0054] In some embodiments, the pupil 210 is tracked using a Purkinje method, which uses a 

camera and light source to compute the eye's orientation based on light reflections from the 
anterior surface of the cornea. Because the Purkinje method does not depend on the pupil 
opening and closing concentrically about the eye's optic axis, the Purkinje method is relatively 
accurate. However, the Purkinje method often requires a stringently-controlled lighting 
environment to be able to detect the rear surface reflection off the eye*s lens. For those 
embodiments that utilize the Purkinje method, the pupillometry image also includes a location of 
the first Purkinje image (hereinafter, "the Purkinje location") 230, in addition to images firom 
multiple eye-illumination light emitting diodes (LEDs) 240. While the Purkinje method is 
specifically disclosed, it should be appreciated that other eye-tracking methods can be employed, 
as should be appreciated by those having skill in the art. 

[0055] FIG. 3 is a diagram illustrating an embodiment of a samphng matrix 320 in a 

refractometer adapted to acquire neuro-ocular wavefi'ont data. Again, for proper orientation, 
both the X-axis and the y-axis are shown in FIG. 3.. The sampling matrix 320 superimposes onto 
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the location of the pupil 210 as identified by the various LED images 230, 240. As shown in 
FIG. 3, in some embodiments, the sampUng matrix 320 comprises muUiple samphng elements 
310a, 310b (hereinafter "sampling elements 310") that are arranged to form a grid-like pattem. 
A treatment envelope 330 within the sampling matrix 320 defines the area for treatment. This is 
shown more clearly in FIG. 4, where the treatment envelope 330 and the sampling matrix 320 are 
shown without the sampling elements 310. 

[0056] FIG. 5 is a diagram showing embodiments of various treatment envelopes 330a . . . 330c 
(hereinafter referred to as "pupil-sampling maps 330"). The left-most image of FIG. 5 shows a 
reduced sampling map 510 that includes a square matrix that is wholly located within the 
treatment envelope. The center image of FIG. 5 shows a nominal sampling map 520, which 
includes sampling points that are substantially within the treatment envelope. The right-most 
image of FIG. 5 shows an increased sampling map 530, which includes any sampling point that 
intersects the treatment envelope. In that regard, it should be appreciated that a sampling matrix 
associated with the pupillometric measurements can be defined in various ways. 

[0057] FIG. 6 is a diagram showing an example video display of a refiractometer adapted to 
acquire neuro-ocular wavefront data. As shown in FIG. 6, the video display includes an eye 
display 610, which shows the pupil tracking box and the Purkinje location with reference to the 
x-y plane. Additionally, the video display comprises a sampling map 620, which shows the 
sampling points for acquiring the neuro-ocular wavefiront data. Additionally, the video display 
includes an instrument control panel 630, which pemiits an operator to adjust various parameters 
associated with pupillometry. For example, for embodiments that employ the Purkinje method, 
the instrument control panel 630 can include a graphical interface for adjusting the Purkinje 
illumination. In some embodiments, the video display can fiirther include an instrument status 
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panel 640, which apprises the operator of whether or not the instrument is properly operating. 
Additionally, the video display can include a tracking control display 660, which can include 
various settings and measurements 650, such as, for example, the measured pupil diameter, the 
pupil position, and the Purkinje offset. While specific display parameters are shown in FIG. 6, it 
should be appreciated that the video display can be configured to include other parameters of 
interest to the operator. Examples of other display parameters are shown in FIG. 7. 
[0058] FIG. 8 is a diagram showing a video display with an alignment target and a test spot. The 
diagram of FIG. 8 is substantially identical to the diagram of FIG. 6, with the exception of the 
tracking display 810. Hence, only the tracking display 810 is discussed in detail with reference 
to FIG. 8. As shown in FIG. 8, the tracking display 810 comprises an alignment target 820 and a 
test spot 830. In the embodiment shown, the alignment target 820 is configured as cross-hairs 
while the test spot 830 is configured as a spot. Specifically, the tracking display 810 represents 
an image that is displayed to the subject. In that regard, for each of the sampling points in the 
sampling map 620, the subject is requested to interactively move the test spot 830 to the center of 
the aUgnment target 820. This process is repeated for each of the sampling points in the 
sampling map 620. Thus, for example, given the 21 -point sampling matrix of FIG. 8, the 
refi-actometer interactively receives the first alignment data fi^om the subject. Upon receiving the 
first alignment data, the next sampling point is selected, and the subject repeats the process for 
the second alignment data. The acquisition of all alignment data, in effect, represents the 
interactive acquisition of the neuro-ocular wavefiront data. From the obtained data, a correction 
is then calculated. In some embodiments, the calculation of the correction can be seen as the 
estimation of one or more correction factors. Once data has been collected fi-om all of the 
sampling points, the collected data can be used to reconstruct the neuro-ocular wavefiront as 
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follows. 

[0059] The neuro-ocular wavefront can be represented as a Zemike polynomial according to: 

Nt 

[0060] Wix,y) = r„,,„£c„ • F„(p,e) [Eq. 1] 

where r„^^ is the normative radius, {p,6) are the Zemike function terms, and C„ are the 

coefficients for the neuro-ocular wavefront and Nt is the number of coefficients used to define W. 
Eq. 1 can be re-written as: 

Nt 

[0061] W{x,y) = Y^B„ F„{p,Q) [Eq. 2] 

where: 

[0062] ^„=r„.™-Q [Eq.3] 

[0063] Given the coordinate system, as defined in FIG. 1, a given sampling point {x,y) on the . 

sampHng map can be defined according to: 
[0064] jc = pcos(e), y = psin(e) [Eq. 4] 

with: 

I 2 2~ 

[0065] p = ifLjtZ_, 0<^x'+y'<r^ [Eq. 5] 

where is the largest radius within the sampling map. 

[0066] As noted with reference to FIG. 1, the degree of aberration at any sampling point can be 
defined by the compensation angles: 

[0067] dx{x,y)^^^ [Eq.6] 

dx 

and: 
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[0068] dy{x,y)=+^ [Eq. 7]. 

dy 

[0069] Since the neuro-ocular wavefiront data can be represented in accordance with Eq. 1, the 

solution of Eq. 1 can be used to determine various neuro-ocular characteristics as well as various 
treatments for neuro-ocular aberrations. In other words, the neuro-ocular system can be 
characterized by solving for the various coefficients in Eq. 1, and these solutions can be used to 
design various treatments. 

[0070] The specific neuro-ocular wavefront of Eq. 1 can be written generically as: 

[0071] f{x)=f^a,X,{x) [Eq.8] 

where X^^ (x) is the basis function for /, and is the coefficient associated with its respective 

basis fimction. Given Eq. 8, the solutions for coefficients can be obtained using a least-squares 
algorithm (or any other fitting routine) such that the merit fimction: 



M 

[0072] x'=2 



[Eq. 9] 



1=1 L *=i 
is minimized. Defining 

[0073] a^dx [Eq. 10] 

and 

[0074] ^^dy [Eq.ll], 

the compensation angles for each sampling point can be represented as: 

Nt Nt 

[0075] a{x,y) = = I.C„G„ [Eq. 12] 

and: 



Page 14 



TKHR 82001-1080 



[0076] 



Nl 



Nt 



[Eq. 13]. 



11=1 



where: 



[0077] 



and: 



[0078] 



n nofTtt 



cbc 



5v 



[Eq. 14] 



[Eq. 15] 



[0079] While Eq. 9 provides the merit function for a single-variable equation, it is clear from 
Eqs. 12 and 13 that the neuro-ocular wavefront is a function of two variables, namely, the 
compensation angle associated with the x-axis and the compensation angle associated with the y- 
axis. Thus, modifying Eq. 9 to accommodate two variables results in: 



[0080] 



Np 



Nt ^ Nt 



-l2 



F,>',) 



[Eq. 16]. 



Solutions to C„ are determined by: 



[0081] 



Np 

1=1 



Nt 
7=1 



7=1 



[Eq. 17] 



where Np is the number of sample points and 
[0082] n€{l,2,:-,N,) 



[Eq. 18]. 



[0083] 



Eq. 17 can be re-written as: 

Np Nl 

i=\ y=i 



Np Nt 

' ^ ' 1=1 j=l 

Np 



[0084] 



1 

i=l 



[Eq. 19] 
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or 



[0085] 



Nt 

z 



( Np 



Nt 



Np 

E 

1=1 



[Eq. 20] 

Simplifying Eq. 20 into a single matrix operation results in: 



[0086] 



Nt 



YS„.Ci=T„ 

y-1 



where: 



[0087] 



Np Np 
i=l /=I 



and: 



[0088] 



Np Np 

1=1 1=1 



[Eq.21] 



[Eq. 22] 



[Eq. 23] 



Using matrix notation, Eq. 22 becomes: 
[0089] S = G^G + H^H [Eq. 24] 

and Eq. 23 becomes: 

[0090] T = G'^a + H^p [Eq. 25] 

[0091] Thus, C„ is derived as a solution to: 

[0092] T = SC [Eq.26] 

where C is a vector that contains the Zemike function terms. Alternatively, Eq. 26 can be 
depicted as: 
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[0093] 



^2.1 ^2,2 



Nt 



. c , 



[Eq. 27] 



to which the solution is: 



[0094] 



C = /«v(S)T 



[Eq. 28] 



or, written differently: 



Nt 



[0095] 



[Eq. 29] 



where: 



[0096] ft = mv(S) [Eq. 30]. 

[0097] Typically, zero-reference measurements are made in order to remove any misalignment 

between the subject and the measurement instrument. The zero-reference is measured by 
calculating a baseline angular deviation associated with a central point at the entrance pupil of 
the subject. The reference measurements can be acquired multiple times and averaged in order 
to increase the accuracy of the zero-reference. The zero-references for the x-axis is given as: 

[0098] = - dx^ [Eq. 31] 

while the zero-reference for the y-axis is given as: 

[0099] Jy,,,,^ = - dy^^^ [Eq. 32] 

[00100] For some instruments, such as the InterWave™ scanner, the neuro-ocular wavefront data 
is acquired with respect to the instrument's sphere-equivalent setting. Thus, for those 
instruments, it is desirable to correct the Zemike coefficients to account for the instrument 
setting. Given the Zemike function terms, as shown in FIG. 31, Eq. 1 can be expanded up to five 
terms as: 
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[00101] 



W(r,Q) = Q •2pcose + C2 '^-o™ •2psine 

+ C3-r„„^ V6 .p^ 008(29) 

+ C,-/-_-V3.(2p^-l) 

+ Q-^„o™-V6-p^ 8111(29) 
•¥ higher order terms 



[Eq. 33]. 



[00102] Since the Zemike function terms associated with the sphere and cylinder are C3, C4, and 

C5, each of these terms can be corrected to account for the instrument's sphere-equivalent setting, 
such that: 



[00103] 



[00104] 



D , r •cos(2-9 .) 

■^cyl 'norm ^'^'^V'^ ^ axt ^ 

4-76 



x-^' ^ ^sph ^norm ^ cyl ^norm 



[Eq. 34] 



[Eq. 35] 



and: 



[00105] 



Dcy,-r„^ •sin(2 -9„,.) 
4-76 



[Eq. 36] 



where D^^^ , D^, , and represent the ophthalmic sphere value, the ophthalmic cylinder value. 



and the ophthalmic axi8 value, respectively, with: 



[00106] 



D 



D 



sph.MR 



sph,CP 



l-TD. 



[Eq. 37] 



sph,m 



and: 



[00107] 



D 



D. 



cyl,m 



eyl,CP 



\-TD 



[Eq. 38] 



where the subscript CP denotes the corneal plane and the subscript MR denotes the eyeglass 
plane. 

[001 08] Given these solutions to Eq. 1 , and their various modifications, an example of a 
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correction obtained using these equations is shown with reference to FIGS. 9 and 10. 
[00109] FIG. 9 is a diagram showing a coordinate system of the neuro-ocular wavefront data and 
a coordinate system for a laser ablation map, to be used to direct laser treatment in vision 
correction that corresponds to the neuro-ocular wavefront data. Specifically, FIG. 9 shows a 
wavefront coordinate system (left) and a coordinate system for laser ablation map (right). The 
wavefront coordinate system and the iablation map coordinate system have a congruent 
relationship to each other. For example, the center of the wavefront coordinate system is the 
pupil centroid 910. On the other hand, the center of the ablation map coordinate system is the 
comeal vertex 920. Thus, if the position of the comeal vertex 920 is known with reference to the 
location of the pupil centroid 910, then an appropriate adjustment can be calculated for the 
ablation map. 

[001 10] For example, presuming that the Purkinje location is coincident with the comeal vertex, if 

the Purkinje offset values are defined according to: 
[001 1 1 ] = PurkinjeOffsetMMx [Eq. 39] 

and: 

[001 12] = PurkinjeOffsetMMy [Eq. 40], 

then the coordinate systems of FIG. 9 are related to each other according to: 

[00113] jc = ^ + jCo [Eq.41] 

and: 

[00114] y^Y^y, [Eq.42]. 

[001 15] Given the relationship between the two coordinate systems, appropriate corrections can 

be calculated. Example procedures for calculating corrections are provided below, with 
reference to the flowcharts of FIGS. 17 through 27B. For simplicity, FIG. 29 provides a 



Page 19 



TKHR 82001-1080 

summary of the variables associated with the calculations in Eqs. 39 through 42. 

[001 16] FIG. 10 is a diagram showing an optical zone 1010 and a transition zone 1020 associated 

with the ablation map of FIG. 9. As shown in FIG. 10, a treatment profile can be represented as 
an aggregate of an optical zone treatment and a transition zone treatment according to: 

[00117] 7'.-r(^,J^)= J^) [Eq. 43]. 

where J^.^^ ^) represents the treatment profile of the optical zone, and T^_j^ (AT, Y) 
represents the treatment profile of the transition zone. Since the first two coefficients, C\ and 
Ca, of Eq. 1 represent tilt of the wavefiront error, which affects the fixation direction of each eye, 
these terms are effectively removed fi-om ablation calculations by setting: 

[00118] C, =0 [Eq.44] 

and; 

[00119] ^=0 [Eq.45]. 

[00120] For some embodiments, it may be desirable to modify the Zemike function terms 

associated with sphere and cylinder in order to apply a different sphero-cylinder correction than 
that predicted by the Zemike function terms themselves. For example, a particular subject may 
indicate a preference for different visual corrections than what would be predicted solely from 
the Zemike function temis. For those cases, the subject's input can be accommodated by 
modifying the C3, C4, and C5 terms, as defined in Eqs. 34 through 36 and, also, as shown in FIG. 
30. Since those modifications are described above, further discussion of such modifications is 
omitted here. 

[00121] Given the definitions of FIG. 30, the treatment profile for the optical zone can be 
represented according to: 
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[00122] [Eq. 46] 

where Rq^ represents the radius of the optical zone, W^.^ represents the minimum value of the 
wavefront error (or the most negative value of the wavefront error) over the area of the optical 
zone, and represents the index of refraction of the cornea. Additionally, from the definitions 
of FIG. 30, the treatment profile for the transition zone can be represented according to: 

[00123] T,_^iXJ) = te +T,.o,{Xo^Joz)) {^- (^''^''''^ X fM i^oz <r^Rn) 

[Eq. 47] 

where represents a constant ablation depth, which is added to faciUtate the addition of the 
transition zone, {XQ^yYoz ) the optical zone treatment amount at a point where the optical 
zone boundary intersects the vector from the pupil center to the point {X, y),.and r is the 
distance from the mesopic pupil center to the point {X, Y) . In Eq. 47, the function /{x^y) is a 
smoothing function that is configured to contour the boundary between the optical zone and the 
transition zone such that there is a smooth transition between these two zones. As such, in some 
embodiments, the smoothing function may be a supergaussian having the form: 

[00124] f(x,y) = e ^''^ [Eq. 48] 

where /i is a non-zero, positive integer value. 
[00125] In some embodiments, Eq. 43 can be modified to include any number of variables, such 
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as, for example, optical parameters, subject parameters, and/or enviromnental parameters, as 
long as those variables affect the treatment in a relatively linear manner. Defining the treatment 
to include such variables as T^.^ {X^ X), that treatment can be defined according to: 

-^ZP.iXJ) [Eq.49] 

k 

where W„ represents the magnitude of the weighting term associated with the variable ^„ , the 
symbol r\ represents material properties, f(^(xj) represents an ablation efficiency compensation 
function, and represents other spatially-dependent, predictable processes. Here, ^„ 

k 

can represent various factors, such as, for example, photopic pupil diameter, mesopic pupil 
diameter, cycloplegic pupil diameter, near-vision preoperative refi"action sphere, near-vision 
preoperative refi-action cylinder, near-vision preoperative refraction axis, far- vision preoperative 
refi-action sphere, far-vision preoperative refi-action cylinder, far-vision preoperative refraction 
axis, near-vision postoperative refi-action sphere, near-vision postoperative refiraction cylinder, 
near-vision postoperative refi-action axis, far-vision postoperative refi-action sphere, far-vision 
postoperative refraction cylinder, far-vision postoperative refraction axis, left eye, right eye, 
asphericity, axis angle, optical zone diameter, transition zone diameter, central pachymetry, 
corneal topographic measurements, spherical aberration as a percent of total root-mean-square 
(RMS) aberration, coma as a percent of total RMS aberration, trefoil as a percent of total RMS 
aberration, high-order aberrations as a percent of total RMS aberration, astigmatism index, 
corneal width, front surface corneal curvature, back surface corneal curvature, front-to-back 
alignment, age, side of dominant eye, preference between day vision and night vision, treatment 
purpose, ethnicity, iris color, gender, temperature, humidity, microkeratome used for comeal 
resection, flap size, time elapsed from opening of flap to ablation, surgeon, estimated total time 
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during opening of flap, expected flap thickness, procedure type, scanner used, laser used, day of 
surgery, location of flap hinge, or any combination of these variables, which are also referred to 
herein as confounding parameters. The coefficients of Eq. 49 can be determined by known 
statistical methods, so long as neuro-ocular wavefront data from a sufficient population of 
subjects has been gathered. Such data can be obtained, for some embodiments, in accordance 
with one or more of the processes set forth in FIGS. 17 through 27B, below. 
[00127] For some embodiments, the ablation efficiency fimction accounts for the loss of 

ablation efficiency across the cornea surface. For those embodiments, the reflectance R^^^ of the 
corneal surface is calculated as: 

[00128] i?^^ = M_LLL [Eq.50] 



where: 



2 



_ «2 cos9, -n, cosBj 



[00129] = ^ -i— J [Eq. 5 1] 

«2 COS9, +«j COS02 

and: 

[00130] ^ /». cose. cose, ^^^ 32] 

W| COS 6 J +/22 COS 9 2 

with «2 representing a complex index of rejfraction for corneal tissue, w, representing a complex 
index of refraction for incident medivun, ^, representing an incident angle, and $2 representing 
an internal angle fefi-acted according to Snell's law. The incident angle is calculated from 
comeal topographic data, Z, at point (X,Y) in accordance with: 

[00131] cos^, "1=T 

thereby resulting in: 
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[00132] 



tane.=J|f 



1 ^ 



[Eq. 54] 



where R is the distance from the corneal vertex to the point of interest (calculated as 



[00133] 
[00134] 



), and (p is the azimuthal angle such that tan ^ = — . 

If the cornea is represented by a conic section of: 
X^+Y^ +Z\l + K)-2ZRo=0 
where AT is a conic constant, and Rq is a paraxial radius of curvature, then: 



[Eq. 55] 



[00135] 



tanO. = 



R_ 
R. 



i-(i.^).|; 



[Eq. 56] 



[001 36] From Eq. 54 or Eq. 56, in conjunction with Eqs. 5 1 through 52, the ablation depth can be 
calculated according to: 



[00137] 



d = do- log, 



1 



[Eq. 57] 



where do is the ablation depth when ^ = e • or ^ = 2.7 1 8 • , <|) is the applied fluence, and 

is the ablation threshold. Thus, for a small beam ablation where the laser bean width is much 
smaller than the pupil diameter, Eq. 57 can be rewritten as: 



[00138] 



c? = rfolog. 



X -(i-^^^J-cose, 



[Eq. 58] 



Given this, a rudimentary approximation of the ablation efficiency can be represented according 
to: 
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[00139] 



log. 



X •(l-^„„J cos0, 



[Eq. 59]. 



If treatment is applied in AT layers, then the transition zone treatment can be represented as: 
[00140] T,_r{X,Y)= Xhi^J) [Eq. 60] 



1 



where Nl(X,Y) represents the number oflayers at point (X.J^.and: 



[00141] 



t(xrt-t d(X,Y) 

h(^>n-tu.yer ^^^ ^^ - hayer 



- (l-R. )-cose., 
9o 



log. 



[Eq. 61] 



with t^^ being a nominal value for the removal rate and ti(X,Y) being the actual material 

removed per layer at that point. If treatment is presumed to be constant across the pupil plane, 
then: 



[00142] 



[Eq. 62] 



layer 



[00143] 



Thus, according to Eqs. 60 and 62: 



[Eq. 63]. 



[00144] FIG. 1 1 is a diagram showing an embodiment of a sampling map that has a predefined 

matrix of sampling points. SpecijScally, FIG. 11 shows a nominal sampling map 520 in which 
the pupil has been divided into a predefined matrix of sampling points. Each sampling point 
represents a visual test zone. The sampling points are shown in FIG. 11 as the small circles of 
various colors. The embodiment of FIG. 1 1 shows a sampling fi-om a subject at a particular 
stage in the data-acquisition process. Specifically, the blank circles 1110 represent those visual 
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zones that have not yet been sampled. The yellow circles 1 140 represent those visual zones that 
have been temporarily skipped and which may be sampled at a later time. These temporarily- 
skipped zones can be skipped as a result of the operator choosing to conduct the data acquisition 
out of order, or, altematively, these zones can be temporarily skipped due to feedback from the 
subject that indicates difficulty in providing input. The dark blue circles 1 130 represent those 
visual zones for which the subject has provided input. The red circles 1 120 represent those 
visual zones that have been permanently skipped due to difficulty in obtaining a reliable reading 
from the subject. The light blue circle 1 150 represents the visual zone for which data is currently 
being obtained. As one can imagine, neuro-ocular wavefront data can be estimated by obtaining 
input from the subject for each of the visual zones. 

[00145] FIG. 12 is a diagram showing an embodiment of a sampling map that permits 

customization of sampling points. For some embodiments, it may be desirable to obtain data 
from various visual zones that do not easily coincide with the predefined sampling map of FIG. 
11. For those embodiments, a customized sampling map 1220 can be used in which predefined 
sampled visual zones 1 130 can be supplemented with custom sample points 1210. By permitting 
customization of the sampling map, more accurate neuro-ocular waveform data can be obtained. 

[00146] FIG. 13 is a diagram showing an embodiment of a subject data entry screen for a 

refractometer adapted to acquire neuro-ocular wavefront data. As described above, in some 
embodiments, subject information can be used to statistically compute various subject-related or 
environment-related factors that may affect treatment- FIG. 13 shows an embodiment of a data 
entry screen by which the subject data can be provided. As shovm in FIG. 13, the data entry 
screen comprises an input 1305 for the gender of the subject, an input 1310 for the date of birth 
of the subject, an input 1315 for the type of examination, and inputs 1320a, 1320b (collectively 
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referred to herein as 1320) for optical characteristics for each eye. The optical characteristics, in 
some embodiments, can include the sphere 1330a, 1330b, the cyHnder 1340a, 1340b, and the 
axis 1 350a, 1350b for each eye. It should be appreciated that the data entry screen can be 
customized, as a matter of design, to include other fields and, also, to omit certain fields. 

[00147] FIG. 14 is a diagram showing an embodiment of a video display that facilitates on-screen 

measurements of physical characteristics of the eye. Specifically, FIG. 14 shows, with reference 
to the X-axis and the y-axis, the pupil tracking box 220 and the Purkinje location 230 as well as 
the eye illumination LED locations 240. The embodiment of FIG. 14 includes an on-screen 
measurement tool 1410 that can be adjusted with a cursor (not shown) using a graphical user 
interface (GUI, not shown). The on-screen measurement tool 1410, for some embodiments, can 
measure the width of the comea. It should be appreciated that the system can be configured to 
accommodate other similar measurements. 

[00148] FIG. 15 is a diagram showing an embodiment of a video display that permits annotation 
of various features of the eye. For this embodiment, various features 1510, 1530 can be selected 
with a cursor (not shown) using a GUI (not shown). The selected features 1510, 1530 can be 
annotated 1520, 1540 in order to provide more information on the characteristics of the eye. 
This information can subsequently be used to refine the diagnosis and treatment of vision. 

[00149] As one can see fi-om FIGS. 14 and 15, the aberrometer that determines the neuro-ocular 

wavefi-ont error can be configured in a variety of ways, thereby permitting greater accuracy in 
diagnosing and treating the visual system of a subject. 

[00150] FIG. 16 is a block diagram showing an embodiment of a system configured to aggregate 

data fi-om multiple spatially resolved refi-actometers. As noted with reference to Eq. 49, 
statistical analysis can be performed on neuro-ocular wavefi*ont data. The embodiment of FIG. 
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16 provides one example, among others, of a system that is configured to obtain multiple sets of 
neuro-ocular wavefront data. 

[00151] As shown in FIG. 16, in some embodiments, multiple clinical sites 1605, 1620 are 
communicatively coupled to a server 1640 over a network 1635, such as, for example, the 
Internet. FIG. 16 also shows a client 1655 that is communicatively coupled to the Intemet 1635. 
Each clinical site 1605, 1620 can include multiple refractometers 1610a . 161 On (collectively 
1610), 1625a . . . 1625n (collectively 1625) that are communicatively coupled to the Intemet 
1635 through routers 1615, 1630 that are located at their respective clinical sites 1605, 1620. 
Since networked environments are known in the art, further discussion of the routers 1615, 1630 
and the Intemet 1635 is omitted here. 

[00152] The server 1640 is also communicatively coupled to the Intemet 1635, thereby permitting 

data exchange between the server 1640 and the various refractometers 1610, 1625. The server 
1640 comprises an information storage device 1645, such as, for example, a high-capacity hard 
drive or other non-volatile memory devices, which are known in the art. Coupled to the 
information storage device 1645 is a processor 1650. The processor 1650, for some 
embodiments, is configured to handle the incoming and outgoing data from the server 1640. 
Since such processor functions are known in the art, further discussion of those functions is 
omitted here. In addition to handling such known functions, the processor 1650 is also 
configured to respond to requests from the client 1655. Thus, for example, the server processor 
1650 receives queries fi-om the client processor 1660 and provides any requested service to the 
cUent 1655 in response to the query. 

[00153] The client processor 1660 is configured to generate requests to the server 1640 and 
receive data in response to those queries. The client processor 1660 is further configured to 
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perform statistical analysis on the requested data. Since those statistical processes, in addition to 
the hardware and software for performing the statistical analysis, are known in the art, further 
discussion of those processor functions is omitted here. 

[00154] As shown in FIG. 16, both the server processor 1650 and the client processor 1660 can be 

configured, using various known software and hardware, to gather data, exchange data, and 
process data. Specifically, for some embodiments, the client processor can be appropriately 
configured to perform the mathematical operations as defined in Eqs. 1 through 63. In other 
words, the processors 1650, 1660 can be configured to generate the neuro-ocular wavefi"ont error 
data using the inputs received from a subject. Additionally, the processors 1650, 1660 can be 
configured to calculate treatment profiles fi"om the neuro-ocular wavefront data. For some 
embodiments, the processor 1650 includes various logic components that are adapted to execute 
the processes described with reference to FIGS. 17 through 27B. Those processes can be 
embodied in a neuro-ocular wavefi-ont program, which computes the neuro-ocular wavefront 
data from the inputs provided by the subject, and statistical regression program, which computes 
the effect of various vision parameters on the neuro-ocular wavefront data. 

[00155] As one can see, by having such a distributed system in which data is stored at a central 
repository 1645, multiple sets of neuro-ocular wavefront data can be aggregated to better 
estimate various effects, such as, for example, subject effects, environmental effects, and optical 
effects. 

[001 56] Having described various embodiments of systems associated with the acquisition and 
analysis of neuro-ocular wavefront data, attention is turned to FIGS. 17 through 27B, which 
show various processes associated with neuro-ocular wavefront data acquisition and analysis. It 
should be appreciated that the following processes, in some embodiments, can be executed by 
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the above-described systems. Alternatively, the processes of FIGS. 17 through 27B can also be 
executed with other devices that are configured appropriately. 

[00157] In a broad sense, the processes of FIGS. 17 through 27B can be seen as an interactive 

process in which a subject (or patient) provides neuro-ocular wavefront data by interacting with a 
spatially-resolved refractometer. For example, in some embodiments, the subject is placed in 
front of a refractometer, and an eye-tracking algorithm tracks the location of the subject's pupil in 
order to eliminate much of the movement artifacts that may be present due to natural eye 
movements. The eye, and more specifically the pupil, is divided into multiple zones. For some 
embodiments, those zones are defined by a matrix of sampling points, which can be predefined 
sampling points or customized sampling points. The refractive characteristic of each zone is 
obtained interactively from the subject. Since that information is obtained with feedback from 
the subject, the obtained data includes neuro-ocular information, rather than merely ocular 
information. Once those characteristics of all of the zones have been obtained, a neuro-ocular 
wavefront can be estimated from those discreet sampling points. As noted above, due to the 
input by the subject, the neuro-ocular wavefront data provides some information on the neuro- 
ocular pathways, rather than merely providing information from the eye alone. 

[00158] FIG. 17 is a flowchart showing an embodiment of a method for acquiring and processing 

neuro-ociilar wavefront data. As shown in FIG. 17, some embodiments of the process can be 
broadly seen as comprising two steps. First, a data acquisition step (1705) and a data analysis 
step (1725). For some embodiments, the data acquisition step (1705) comprises the step of 
obtaining (1720) neuro-ocular wavefront data from a subject, while the data analysis step (1725) 
comprises the step of ascertaining characteristics of a visual system from the obtained neuro- 
ocular wavefront data. The neuro-ocular wavefront data can be obtained (1720), in some 
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embodiments, using one or more of the systems described with reference to FIGS. 1 through 16. 

[00159] In other embodiments, the data analysis step (1725) comprises the step of correlating 
(1730) neuro-ocular wavefront data to parameters associated with the visual system of the 
subject. The parameter associated with the visual system can be, but is not limited to, an optical 
parameter, a subject parameter, or an environment parameter. The optical parameter, the subject 
parameter, and the environmental parameter, in the aggregate, are referred to herein as "vision 
parameters." In some embodiments, the optical parameters can include, but are not limited to, 
photopic pupil diameter, mesopic pupil diameter, cycloplegic pupil diameter, near-vision 
preoperative refraction sphere, near- vision preoperative refraction cylinder, near-vision 
preoperative refraction axis, far-vision preoperative refraction sphere, far-vision preoperative 
refraction cylinder, far-vision preoperative refraction axis, near-vision postoperative refraction 
sphere, near-vision postoperative refraction cylinder, near- vision postoperative refraction axis, 
far-vision postoperative refraction sphere, far-vision postoperative refraction cylinder, far-vision 
postoperative refraction axis, left eye, right eye, asphericity, axis angle, optical zone diameter, 
transition zone diameter, central pachymetry, spherical aberration as a percent of total root- 
mean-square (RMS) aberration, coma as a percent of total RMS aberration, trefoil as a percent of 
total RMS aberration, high-order aberrations as a percent of total RMS aberration, astigmatism 
index, comeal width, front surface comeal curvature, back surface comeal curvature, front-to- 
back alignment, or any combination of these optical parameters. 

[00160] The subject parameters, in some embodiments, can include, but are iiot limited to, age, 
side of dominant eye, preference between day vision and night vision, treatment purpose, 
ethnicity, iris color, gender, or any combination of these subject parameters. 

[00161] For some embodiments, the environmental parameters can include, but are not limited to. 
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temperature, humidity, microkeratome used for corneal resection, flap size, time elapsed from 
opening of flap to ablation, surgeon, estimated total time during opening of flap, expected flap 
thickness, procedure type, scanner used, laser used, day of surgery, location of flap hinge, or any 
combination of these environmental parameters. 

[00162] Thereafter, the data analysis step (1725) proceeds to calculate (1740) a correction 

corresponding to one or more treatments for reducing the anomalies in a visual system of the 
subject. The calculation of the correction can be seen as a calculation of one or more correction 
factors. Thus, the correction can correspond to a prescription for spectacles. In other 
embodiments, the correction can correspond to a prescription for contact lenses. However, in a 
preferred embodiment, the correction can correspond to a treatment for a refractive surgical 
technique, such as, but not limited to, radial keratotomy (RK), astigmatic keratotomy (AK), 
automated lamellar keratoplasty (ALK), photorefractive keratectomy (PRK), laser in situ 
keratomileusis (LASIK), intracomeal ring segments (Intacs), intracomeal lens insertion, laser 
thermal keratoplasty (LTK), phakic intraocular lenses, or any coinbination of these refractive 
surgical techniques. Likewise, the correction factors can correspond to one or more factors that 
affect an ablation profile for a corresponding refractive surgical method. 

[00163] FIG. 18 is a flowchart showing, in greater detail, an embodiment of the data-acquisition 

step (1720) from FIG. 17. As shown in FIG. 18, the step of obtaining (1720) the neuro-ocular 
wavefront data: can be seen, for some embodiments, as a two-step process. In the first step, the 
process identifies (1820) visual zones corresponding to different regions of the eye. Thereafter, 
the process interactively obtains (1830) information related to the visual zones. The information 
can be obtained using an interactive refractometer, such as, for example, one or more 
refractometers that have been described with reference to FIGS. 1 through 16. 
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[00164] FIG. 19 is a flowchart showing, in greater detail, an embodiment of the step of 

identifying (1820) visual zones, from FIG. 18. As shown in FIG. 19, in some embodiments, the 
process identifies (1920) an area associated with the entrance pupil of the eye. Thereafter, the 
process overlays (1930) a virtual matrix onto the identified area. Preferably, the area associated 
with the entrance pupil of the eye is identified and tracked using the Purkinje method, as 
described above. Thus, when the virtual matrix is overlayed (1930) onto the pupil area, that 
matrix will be in substantially the same place, regardless of spurious eye movements. 

[00165] FIGS. 20A and 20B are flowcharts showing, in greater detail, an embodiment of the step 
of interactively obtaining (1830) information, from FIG. 18. As shown in FIG. 20A, an 
embodiment of the process begins (2005) when a reticule image is projected (2010) at the center 
of the eye and data is gathered from the various sampling points or visual zones. In that regard, 
the process determines (2015) whether or not all visual zones have been analyzed and, so long as 
all of the visual zones have not been analyzed, further determines (2025) whether or not a 
customized sampling map should be used for the analysis. If the process determines (2025) that 
the analysis is customized, then the next customized visual zone is selected (2030). For some 
embodiments, this can be done interactively with the operator of the refractometer through a 
graphical user interface (GUI) similar to that shown in FIG. 12. If, on the other hand, the process 
determines (2025) that the analysis is not customized, then the next visual zone is selected (2035) 
from a predefined matrix. In some embodiments, the predefined matrix can be configured 
similar to the matrix shown in FIG. 11. Regardless of whether a predefined visual zone is 
selected (2035) or a customized visual zone is selected (2030), the process continues to FIG. 
20B. 

[00166] As shown in FIG. 20B, a target image is projected (2050) at the selected visual zone. 
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Upon projecting the target image, the subject is queried (2055) for input, and that input is 
received (2060) from the subject. In some embodiments, the query can be configured as a 
display having an alignment target and a test spot, similar to the display shown in FIG. 8. The 
input from the subject can be received, for some embodiments, using a joystick or other input 
mechanism that permits the subject to aHgn the test spot with the alignment target. Upon 
receiving (2060) the input from the subject, the received input is then stored (2065), and the 
process retums to FIG. 20A to determine (2015) whether or not all visual zones have been 
analyzed. In the event that all visual zones have been analyzed, the process exits (2020). 
[00167] A blur image (also referred to herein as a blur anatomy) is one possible representation of 

image quality of human vision, insofar as the blur image has the following characteristics. First, 
the blur image plainly describes something that is representative of the quality of a subject's 
vision, and it is something that the subject can relate to in simple terms. Second, it provides 
some figure of merit that can be optimized using the parameters that are available to a surgeon as 
control variables. Also, there is a high degree of correlation between the simulation of the 
optimized condition and the actual experience of the patient when treated with the same 
conditions. 

[00168] The blur image depicts a composite image perceived by the patient viewing a point 

source at an infinite distance, such as a bright star or a street light at a great distance. In one 
particular embodiment, the pupil is divided into three circular zones having outer diameters of 
approximately 3.65 mm, 5.16 mm, and 6.32 mm. Rays from an object are traced through the 
pupil and then to the image plane where they are color coded according to the pupillary zone 
through which they pass. Note that the image plane is determined with spatially resolved 
refractometry by localizing the patient preferred retinal plane and taking into account the effects 
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of neural processing. This depiction provides a simultaneous estimate of vision quality for three 
different pupil diameters. The blur anatomy diagnostic is very similar to a conventional "spot 
diagram" used in optical system engineering, with the exception that the color coding provides 
additional information for different pupil diameters. 
[00169] The location where a ray is incident on the retina can be calculated from the wavefront 
error using: 

[00170] a{x,y) = ^— [Eq. 64] 

cbc 

and: 

[00171] ^{^^y)=^^ [Eq.65] 

dy 

where x, and y are the horizontal and vertical coordinates in the entrance pupil plane. In this 
formulation, the coordinate position actually represents the perceived location of the ray in object 
space. Thus, Eqs. 64 and 64 are measured in radians or milliradians. 
[00172] Once the plurality of points in the blur anatomy is calculated using Eqs. 64 and 65, the 
root-mean-square (RMS) blur radius of the blur pattern can be calculated using a standard 
statistical method. The first step in performing this calculation is to determine the coordinates of 
the blur centroid as: 

[00173] oi = -hr- [Eq.66] 

N 



and: 



[00174] P=-V [Eq.67] 

N 
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Once the coordinate of the centroid is determined, the calculation of RMS blur radius is 
completed according to: 



[00175] RMS blur radius = J-^ — [Eq. 68]. 

[001 76] Corrections over the entire pupil are determined in the following manner. First the pupil 
is divided into several annular zones. Next, each of these zones is independently optimized 
using a conventional optimization algorithm. This calculation takes as inputs, the sphere, 
cylinder, and axis of a hypothetically-applied correction and determines the resultant RMS bliw 
radius as an output figure of merit. The numerical optimization proceeds by finding the sphere, 
□, cylinder, □, and axis, that minimizes the RMS blur diameter using a representation 
according to: 

[00177] RMS blur radius = /(9,y,e) [Eq. 69] 

[001 78] These corrections relate to the neuro-ocular wavefi-ont error in accordance with Eqs. 33 
through 36. 

[00179] In this embodiment of optimizing vision, there is variabiUty associated with the method 
used to select the zone diameters. One embodiment of a method used to determine the zone 
diameters is denoted herein as the "equal area" method. This method endeavors to select the 
zone sizes by specifying that the contribution of each zone with relation to the image brightness 
should be equal for each zone. If the pupil diameter under pharmacological dilation is presumed 
to be the maximum pupil diameter possible, then the diameter of the i*^ zone, Di, can be 
determined firom the maximum pupil diameter, Dmax^ and the nxmiber of zones, iV, using: 
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N-2 



[00180] ^"-^ " N [Eq.70] 

etc. 
etc. 



Thus, for example, when A^= 3, the outer diameters of the three zones is: 

[00181] 



D,=D^, D, =^™x • j|. and A =Z)_ [Eq. 71] 

[00182] One consequence of this approach to zone size selection is that the peripheral zones get 
progressively thinner. Besides having the effect of equalizing the intensity contribution of each 
zone, this method is also amenable for optimizing zone diameters of eyes with third-order 
spherical aberration. The reason being that, in this approach, zones designed according to Eq. 
70 will become narrow at the same rate that the variation in sphere error increases. Therefore, 
each zone defined using this method will be treated with substantially the amount of spherical 
correction. 

[00183] Another and more comprehensive method of selecting the zone diameters would be to 
utilize the same optimization technique used to optimize sphere, cylinder, and axis, except that 
zone diameter now becomes one of the variables as an input to the blur diameter calculation. In 
that optimization, a more complicated merit function can be designed, such that the average blur 
for each zone is quantified. This merit function may be designated "total vision blur," and for the 
case of three optical zones can be calculated according to: 
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Total Vision Blur = 



00184] y vvi./P iJ [Eq.721. 

The total vision blur represents the sum of the RMS blur for all three zones. 
[00185] Although the optimization of this function containing twelve variables is considerably more 

complicated than each of the independent zone optimizations, such problems are relatively 

insignificant, given the available processing power. 
[001 86] In physiological optics, it is known that the outer or peripheral parts of the pupil 

contribute less to visual stimuli than do those from the central parts of the pupil. This effect, 

know as the Stiles Crawford effect, typically causes peripheral rays to contribute 

commensurately less to the perception of blur. This factor can be included in the calculation of 

RMS blur radius by modifying Eqs. 66 through 68 according to: 

[00187] OL=^ [Eq.73] 

I 

[00188] P=-^v^ [Eq.74] 

I 

and: 



En.[(a.-a)^+(P,-P)^] 
[00189] RMS blur radius^ M [Eq. 75]. 

[00190] As before, the sum over / is carried out over all the points traced through the pupil from 
the object to the retina. As described above, such processes allows for computation of sphero- 
cylindrical corrections through blur optimization. 
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[00191] FIG. 21 is a flowchart showing, in greater detail, an embodiment of the step of correlating 
(1730) the neuro-ocular wavefront data to vision parameters, from FIG. 17. Specifically, FIG. 
21 shows an embodiment in which the blur image is obtained. The blur image provides an 
emulation of the degree of blur that a subject sees without corrected vision. Additionally, 
depending on the parameters of interest, the blur image can also provide information of how a 
particular subject sees, should that subject's vision be corrected using various treatments. As 
shown in FIG. 21, the process begins by computing (2120) a pupil function from the neuro- 
ocular wavefront data. The pupil fiinction can be calculated in accordance with: 



Pix,y) = T{x,y)-e ^ , V^'+/<- 



pupil 



[00192] 2 



2 

where □ is the wavelength of incident light, Dpupu is the diameter of the entrance pupil, and 
T(x,y) represents an apodizing function that can be used to model variations in effective pupillary 
transmission. For example, the Stiles-Crawford effect can be simulated according to: 

[00193] Tix.y) = e'^'"^'"^ [Eq. 77] 

with a being 0.16 (according to Stiles) or 0.105 (according to Moon and Spencer). 

[001 94] Upon computing (2120) the pupil function, a point-spread function (PSF) is obtained 

(2130) by Fourier transformation of the neuro-ocular wavefront data. The PSF can be computed 
according to: 

[00195] /p5F(a,P) = Co -|F(w,v)| ' [Eq. 78] 

ii=a/X,v=p/X 

where (□, □. represents the angular position of a point in the far-field, F(u,v) is defined as the 
Fourier Transform of P(x,y), and Co is a constant. 
[00196] A color image is resolved (2140) into its component color images (e.g., red, green, and 
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blue images). Each of the component color images is convolved (2150) with the PSF, thereby 
generating component blurred images. The component blurred images are then combined to 
reconstruct (2160) the blur image. For some embodiments, the blur image can be generated by 
the processor 1660 of FIG. 16, so long as the processor is configured with the appropriate 
software and hardware to perform such a function. Since those having skill in the art should 
know how to appropriately configure the processor (e,g., write software code for a particular 
hardware platform), fiirther discussion of such hardware and software is omitted here. 
[00197] FIGS. 22A and 22B are flowcharts showing, in greater detail, an embodiment of the data- 
analysis step (1725) fi-om FIG. 17. As shown in FIG. 22 A, some embodiments of the data- 
analysis step (1725) begin when acquired data is retrieved (2210) and a model equation is 
obtained (221 5). The acquired data represents the neuro-ocular waveform, while the model 
equation represents, for some embodiments, a Zemike polynomial. For other embodiments, the 
model equation is an expansion of the Zemike polynomial, similar to that shown in Eq. 33. 
Upon acquiring (2210) the data and obtaining (2215) the model equation, the tih terms of the 
Zemike polynomial are set (2220) to zero, as shown with reference to Eqs. 44 and 45. 
Thereafter, the sphero-cylinder terms are adjusted (2225) in accordance with Eqs. 34 through 36. 
Once the sphero-cylinder terms have been adjusted (2225), an optical zone ablation profile is 
computed (2230). For some embodiments, this is done in accordance with Eq. 46. The process 
continues to FIG. 22B. 

[00198] As shown in FIG. 22iB, in addition to computing (2230) the optical zone ablation profile, 
the process fiirther computes (2245) a transition zone ablation profile. For some embodiments, 
the transition zone ablation profile is computed (2245) in accordance with Eq. 47. The optical 
zone ablation profile and the transition zone ablation profile are then combined (2250) to 
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produce a total ablation profile. Upon producing the total ablation profile, the process 
determines (2255) whether or not the boundary between the two zones is smooth. If the 
boundary is not smooth, then the transition zone ablation profile is recomputed (2260), and that 
transition zone ablation profile is combined (2250) with the optical zone ablation profile. When 
the boundary between the two zones is smooth, the effect of environmental factors is computed 
(2265). For some embodiments, the effect of those factors can be represented by Eqs. 49 
through 63. The total ablation profile is then adjusted (2270) to correct for those factors. That 
ablation profile is then stored (2275). In some embodiments, the profile is stored (2275) locally 
while, in other embodiments, the profile is stored (2275) at a remote repository. As shown in the 
embodiment of FIGS. 22A and 22B, such a process provides a different approach to vision 
correction. 

[00199] FIG. 23 is a flowchart showing, in greater detail, another embodiment of the data-analysis 
step (1725) firom FIG. 17. As shown in FIG. 23, for some embodiments, the data-analysis step 
(1725) can be seen as a two-step process in which an equation is generated (2310) firom the 
stored inputs, and correction factors are calculated (2315) by inverting the equation. In some 
embodiments, the process of FIG. 23 can be seen as the computation of the coefficients in 
accordance with one or more of Eqs. 1 through 38. 

[00200] FIG. 24 is a flowchart showing another embodiment of a method for acquiring and 
processing neuro-ocular wavefiront data. As shown in FIG. 24, for some embodiments, the 
process begins with psychometric testing being performed (2410) on a subject (or patient). Since 
various psychometric testing is known in the field, fiirther discussion of psychometric testing is 
omitted here. Upon performing (2410) the psychometric tests, psychometric parameters are 
derived (2415). In addition to performing psychometric tests, neuro-ocular wavefront data is 
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obtained (2420) from the subject. That neuro-ocular wavefront data is fit (2425) to an equation. 
For some embodiments, the equation is a Zemike polynomial that represents a neuro-ocular 
wavefront. From the equation, various predicted vision parameters are derived (2430). Upon 
deriving (2430) the predicted vision parameters, the process determines (2435) whether or not 
the vision parameters are similar to the psychometric parameters. If there is a substantial 
discrepancy between the vision parameters and the psychometric parameters, then the fitting 
parameters for the equation are adjusted (2440) and the neuro-ocular wavefront data is fit (2425) 
using the adjusted parameters. If the vision parameters and the psychometric parameters are in 
substantial agreement, then the coefficients of the equation are stored (2445). 

[00201] FIG. 25 is a flowchart showing an embodiment of a method for acquiring and processing 
neuro-ocular wavefront data, in which the processing of subsequently-acquired data is refined. 
As such, the embodiment of FIG. 25 can be seen as a process in which preoperative neuro-ocular 
wavefront data is interactively obtained (2510) from a subject. From the preoperative neuro- 
ocular wavefront data, a correction is determined (2515), and a result of applying the correction 
is predicted (2520). Once the subject has been treated, in accordance with one or more available 
treatments, as outlined above, postoperative neuro-ocular wavefront data is interactively obtained 
(2525) from the subject. Thereafter, the deviation between the preoperative neuro-ocular 
wavefront data and the postoperative neuro-ocular wavefront data is determined (2530). From 
the determined deviation, future corrections are calculated (2535). As seen from the 
embodiment of FIG. 25, by having both postoperative and preoperative neuro-ocular wavefront 
data, future treatments can be refined. 

[00202] FIG. 26 is a flowchart showing an embodiment of a method for acquiring and processing 
neuro-ocular wavefront data, in which multiple sets of data are statistically analyzed. As shown 
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in FIG. 26, multiple sets of neuro-ocular wavefront data are interactively obtained (2620). For 
some embodiments, the multiple sets of data can be obtained from multiple clinical sites, such as 
those shown in FIG. 16. The interactively-obtained neuro-ocular wavefront data is then stored 
(2630). For some embodiments, the data is stored locally, while, for other embodiments, the data 
is stored remotely. Thereafter, the stored neuro-ocular wavefront data is statistically analyzed 
(2640). For some embodiments, the statistical analysis can be seen as a two-step process having 
the steps of applying (2650) a statistical regression across the multiple sets of neuro-ocular 
wavefront data, and determining (2660) parameters of a fitting algorithm using the results of the 
statistical regression. The fitting algorithm can be a least squares fit of the coefficients of a 
Zemike polynomial, as described above. As previously noted, the aggregation of multiple sets of 
data provides an approach to refining treatment methods, insofar as effect from similar 
treatments can be signal averaged over the multiple sets of data. 
[00203] FIGS. 27A and 27B are flowcharts showing another embodiment of a method for 
acquiring and processing neuro-ocular wavefront data, in which multiple sets of data are 
statistically analyzed. In this embodiment, the process begins with the interactive obtaining 
(2710) of preoperative neuro-ocular wavefront data from an initial subject. Thereafter, an initial 
correction is calculated (2715) from the preoperative neuro-ocular wavefront data of the initial 
subject. Also, an initial result of applying that correction is predicted (2720). This can be done 
in accordance with one or more of the above-identified Eqs. 1 through 75. Once the initial 
subject has been treated, in accordance with one or more of the above-identified treatments, 
postoperative neuro-ocular wavefront data is interactively obtained (2725) from the initial 
subject. Subsequently, a deviation between the postoperative neuro-ocular wavefront data and 
the preoperative neuro-ocular wavefront data is determined (2730) for the initial subject. The 
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process continues to FIG. 27B, where the various coefiicients of the corresponding equations are 
adjusted (2745) as a function of the determined deviation. It should be appreciated that the 
coefficients are adjusted as needed. Thus, for some embodiments, only a subset of the 
coefficients are adjusted to account for substantial deviations from the predicted values. 

[00204] The process of FIG. 27B continues by interactively obtaining (2750) preoperative neuro- 
ocular wavefront data for a subsequent subject. Thereafter, a correction is calculated (2755) for 
that subject using the adjusted equation. The process then predicts (2760) a result of applying 
that calculated (2755) correction. Once the subsequent subject has been treated, postoperative 
neuro-ocular wavefront data is interactively obtained (2765) from that subject, and the deviation 
between the predicted result and the postoperative neuro-ocular wavefront data for that subject is 
determined (2770). From that determined (2770) deviation, the process further determines 
(2775) whether or not the adjustment to the equation was appropriate. If the adjustment is 
determined to be appropriate, then the process ends, and the derivation of the coefficients for the 
equation has been sufficiently refined. Alternatively, if the adjustment is determined to be 
insufficient or too extreme, then the process recursively readjusts the coefficients of the equation 
until an appropriate adjustment has been made. 

[00205] As shown with reference to FIGS. 1 through 31, visual aberrometry is an improvement to 

both visual acuity and ocular aberrometry, insofar as visual aberrometry uses neuro-ocular 
wavefront data. As such, visual aberrometry determines the patient preferred retinal plane and 
takes into accoimt the neurological pathways by combining subject feedback with the physio- 
optical characteristics of the eye. The aggregate effect of the neurological pathway, in 
combination with the characteristics of the eye, provides more accurate information for vision 
diagnosis, vision treatment, or a combination of both. 



Page 44 



TKHR 82001-1080 

[00206] The processor, and the components that perform the various functions of the processor, 
may be implemented in hardware, software, firmware, or a combination thereof In the preferred 
embodiment(s), the processor and its functional components are implemented in software or 
firmware that is stored in a memory and that is executed by a suitable instruction execution 
system. If implemented in hardware, as in an alternative embodiment, the processor and its 
functional components can be implemented with any or a combination of the following 
technologies, which are all well known in the art: a discrete logic circuit(s) having logic gates 
for implementing logic functions upon data signals, an application specific integrated circuit 
(ASIC) having appropriate combinational logic gates, a programmable gate array(s) (PGA), a 
field progranmiable gate array (FPGA), etc. 

[00207] Any process descriptions or blocks in flow charts should be understood as representing 
modules, segments, or portions of code which include one or more executable instructions for 
implementing specific logical functions or steps in the process, and altemate implementations are 
included within the scope of the preferred embodiment of the present invention in which 
functions may be executed out of order from that shown or discussed, including substantially 
conciurently or in reverse order, depending on the functionality involved, as would be 
understood by those reasonably skilled in the art of the present invention. 

[00208] The neuro-ocular wavefront program and the statistical regression program, which 

comprise an ordered listing of executable instructions for implementing logical functions, can be 
embodied in any computer-readable medium for use by or in connection with an instruction 
execution system, apparatus, or device, such as a computer-based system, processor-containing 
system, or other system that can fetch the instructions from the instruction execution system, 
apparatus, or device and execute the instructions. In the context of this document, a "computer- 
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readable medium" can be any means that can contain, store, communicate, propagate, or 
transport the program for use by or in connection with the instruction execution system, 
apparatus, or device. The computer-readable medium can be, for example but not limited to, an 
electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, 
device, or propagation medium. More specific examples (a nonexhaustive list) of the computer- 
readable medium would include the following: an electrical connection (electronic) having one 
or more wires, a portable computer diskette (magnetic), a random access memory (RAM) 
(electronic), a read-only memory (ROM) (electronic), an erasable programmable read-only 
memory (EPROM or Flash memory) (electronic), an optical fiber (optical), and a portable 
compact disc read-only memory (CDROM) (optical). Note that the computer-readable medium 
could even be paper or another suitable medium upon which the program is printed, as the 
program can be electronically captured via, for instance, optical scanning of the paper or other 
medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and 
then stored in a computer memory. 
[00209] Although exemplary embodiments have been shown and described, it will be clear to 
those of ordinary skill in the art that a nimiber of changes, modifications, or alterations to the 
invention as described may be made. All such changes, modifications, and alterations should 
therefore be seen as within the scope of the disclosure. 
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